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ABSTRACT 

Context. The Galactic centre is a bright y-ray source with the GeV-TeV band spectrum composed of two distinct 
components in the 1-10 GeV and 1-10 TeV energy ranges. The nature of these two components is not clearly understood. 
Aims. We investigate the y-ray properties of the Galactic centre to clarify the origin of the observed emission. 
Methods. We report imaging, spectral, and timing analysis of data from 74 months of observations of the Galactic 
centre by Fermi/LAT y-ray telescope complemented by sub-MeV data from approximately ten years of INTEGRAL/ 
PICsIT observations. 

Results. We find that the Galactic centre is spatially consistent with the point source in the GeV band. The tightest 3 ct 
upper limit on its radius is 0.13° in the 10 — 300 GeV energy band. The spectrum of the source in the 100 MeV energy 
range does not have a characteristic turnover that would point to the pion decay origin of the signal. Instead, the source 
spectrum is consistent with a model of inverse Gompton scattering by high-energy electrons. In this a model, the GeV 
bump in the spectrum originates from an episode of injection of high-energy particles, which happened ~ 300 years ago. 
This injection episode coincides with the known activity episode of the Galactic centre region, previously identified using 
X-ray observations. The hadronic model of source activity could be still compatible with the data if bremsstrahlung 
emission from high-energy electrons was present in addition to pion decay emission. 
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1. Introduction 

The Galactic centre (GC) is a unique astronomical source 
hosting the nearest supermassive black hole of mass M ~ 
4 X 10®Mq accreting matter in a radiatively inefficient 


Hasselwander et ah] 1 1998 


20071. 


■g/s (Genzel et al. 11994 

Mayer- 

Aharonian et al.| 2004 

|Melia 


The X-ray and infrared source could be associated with 
the supermassive black hole because of its variable nature 
|Baganoff et al. |2001t |Genzel et al.|2003|) . The source of the 
GeV and TeV y-rays is not variable and it is not clear if the 
emission comes from the vicinity of the black hole, if it is 
produced in an extended region around the black hole, or, fi¬ 
nally, if it comes from an unrelated source positionally coin¬ 
cident with the supermassive black hole in Sgr A*. The best 
positional uncertainty of 6” on TeV source HESS J1745-303 
(Acero et al.|2010|) still allows several sources to be respon- 
sible for the observed ver y high energy(VHE) emission: the 
central black hole itself ( Aharonian fc Neronov 2005a|), a 
plerion within several arcseconds from the GG (Wang et al. 
2006 Hinton fc Aharon ian|200 7 ), a putative black hole ple- 
rion around Sgr A’^ ( Atoyan fc Dermer ||2004 ), and the cen¬ 
tral diffuse region around the GG ([Aharonian fc Neronov 


2005b Ballantyne et al.||2007 Ghernyakova et al.||2011[).' 

In the GeV energ y range, the GG sou rce is localised with 
an accuracy of 72” (Nolan et al. 2012) by the Large Area 


Teles cope (LAT) on board the Fermi satellite (Atwood et al. 
2009). An even larger range of possible counterparts could 
be discussed in this energy band. In addition to the as- 
trophysical explanations, the observed GeV emission might 
also originate from the annihilation of dark matter particles 


see e.g. Hooper fc Linden 2011 Hooper & Goodenough 


2011 Abazajian et al. 2014). In this scenario, the signal 


from the central “spike” of the Galactic dark matter halo 
is expected to be extended. 

The analysis of 25 mon t hs of Fermi/LAT observa¬ 
tions by Ghernyakova et al. (2011) has revealed a GeV 
“bump” in the source spectrum. Assuming the identihca- 
tion of the Fermi sourc e with the GG and HESS J1745- 
303, Ghernyakova et al. (2011) proposed a model that al- 
lowed the authors to explain the combined Fermi/LAT - 
HESS spectrum with a model of tt® decay y-ray emission 
from protons diffusing away from the supermassive black 
hole. A similar model was considered by|Fatuzzo fc Melia 


(2012). An alternative explanation of the GeV-TeV spec- 
trum with lepto nic models was proposed by [Kusunose fc:| 
Takahara (2012). Within this model, the GeV component 
of the spectrum is explained by the inverse Gompton (IG) 
emission from electrons injected during NIR/X-ray flares 
of the central source, while the TeV emission originates 
from an unrelated source, e.g. the PWN G359.95-0.04. Still 
another hybrid leptonic-hadronic model ( |Guo et ah 2013) 
combines both approaches, explaining the TeV emission as 
caused by collisions of hadronic cosmic rays (GR) with the 
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Fig. 1. Test-statistic (TS) maps in 60-100 MeV, 100-300 MeV, and 300 MeV-1 GeV energy bands. Only diffuse and 
point-like sources beyond 0.7° from the position of GC were subtracted. The maximum of TS value corresponds to the 
most probable position of added point-like source. The contours for the localisation significance are shown with green 
solid lines for 60-100 MeV map and correspond to significances of 1, 2, and 3cr. The 2FGL catalogue coordinates of the 
GC only deviate from the maximum TS position at ^ Icr level. The position of maximum shifts with the increase of 
energy to the position of Galactic centre. 


surrounding gas and attributing the GeV flux to IC emis¬ 
sion of electrons accelerated near the GC. 

In this paper, we report an updated study of the spec¬ 
tral, timing, and imaging characteristics of the GC source 
in the energy range 60 MeV - 300 GeV, using six years of 
observations with Fermi/LAT. We discuss the implications 
of our findings for the theoretical models. 

2. Data analysis and results 

2.1. Fermi/LAT data analysis 

We use 74 months of Fermi/LAT data (from August, 4th, 
2008 to September, 25th, 2014). We perform the binned 
likelihood analysis, using the most recent Fermi Science 
Tools Software v9r33p0 and the P7REP response functions. 

The analysis is based on the fitting of selected model of 
the observed sky region to the data. We consider a region 
of radius 13° around the GC and P7REP.SOURCE event class 
selectioij^for analysis at energies above 500 MeV. The spa¬ 
tial model includes diffuse Galactic and extragalactic back¬ 
grounds and all sources from the two-year Fermi catalogu^^ 
(2FGL, iNolan et al.||2012|) . 

In the energy band below 500 MeV we use a broader, 
30°, region and higher quality P7REP.CLEAN class pho¬ 
tons. The larger size of the analysis region and a more 
restricted event selection are necessary because of a sig¬ 
nificant increase of the Fermi/LAT PSF (95% containment 
for SOURCE class photons, which are ^ 3° at 500 MeV and 
^ 10° at 100 MeV). To avoid a significant increase in the 

^ http://fermi.gsfc.nasa.gov/ssc/data/analysis/documentation/ 
Cicerone/Cicerone_L AT _IRFs/IRF _overview.html 

^ When the paper was ready for submission, the four-year 
(3FGL) Fermi catalogue was released. In this catalogue, the 
emission is split into two sources located within < 0.1° central 
region. We verified, using the 3FGL catalogue, that our analysis 
is valid for the total emission from this region. 


computational time for fitting the data, as a result of the 
large number of catalogue sources in the selected region, we 
include only sources with a 100 MeV flux higher than 10% 
of the GC. This flux limit is comparable to the systematic 
uncertainties of the flux measurement^ 


2.1.1. Fermi/LAT imaging analysis 


For imaging analysis we build test statistics (TS) maps in 
the energy bands 60-100 MeV, 100-300 MeV, and 300 MeV- 
1 GeV, see Fig. 

To determine the best-fit position of the source, we re¬ 
move all sources beyond 0.7° from the GC position from 
the region model. In the lowest (60-100 MeV) energy band, 
we find that the maximum of TS is shifted by 0.2° from 
the GC. This deviation decreases with energy, and in the 
300 — 1000 MeV band the source position coincides with 
the GC with an accuracy better than the pixel size (0.05°). 
Note, however, that the observed deviation Is not statisti¬ 
cally significant, see Fig. [^left panel. The position of the 
GC lies just beyond the Icr contour. 

At high energies (> 10 GeV), the Fermi/LAT PSF is 
quite narrow (~ 0.2° 68% PSF containment). We start with 
this energy band to put the tightest possible constraints 
on the GC spatial extent. We replace the GC polnt-llke 
source with a uniform disk. We first perform the fitting 
procedure for the model without the GC, and afterwards 
add the disk with variable radius, representing a diffuse GC 
source, to the model. The significance of the detection of 
the disk is given by the difference between log-like values 
of these models (see e.g. Mattox et al.|199^. The fol l owing 

1 inINola ' 


procedure is similar to that described : 


Ian et al. (2012) 


for the localisation of catalogue sources. 

Fig. ^ shows the results of the Imaging analysis. In each 
energy band, the TS value remains roughly the same for 


See http://fermi.gsfc.nasa.gov/ssc/data/analysis/LAT_caveats.html 
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Fig. 2. The change (in the comparison with the point source) of the TS value attributed to the GC modelled as a disk 
of certain radius. The increase of —ATS corresponds to the overall worsening of the fit. The change of —ATS by 9 
corresponds to the Scr upper limit on the size of the disk. The flux, attributed to the disk is shown with red points with 
errorbars. Bottom left panel: The spectra of 1° radius template of the GC, and GC as a point source. Bottom right panel: 
3a upper limit on the size of the disk as a function of energy. The fit with the power law {R = 0AA°{E/\GeV)~'^'^^) is 
shown with green dot-dashed line. 
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Fig. 3. Left panel: Combined Fermi/LAT -HESS spectrum of the Galactic centre with added systematic errors. The 
spectrum of the Galactic bulge from “constant rectangle” model (renormalised to match the GC flux in the highest energy 
Fermi data point) is shown with grey shaded region, see text for the details. Right panel: the flux in 60-100 MeV energy 
bin as a function of minimal catalogue flux of the sources included into the model(for LRYusif ovXCD4z6R30_Tsl50jmag2 
model of the galactic diffuse emission). The corresponding results for gll_iem_v05_revl. f its, gll_iem_v05. f its and 
“constant rectangle” diffuse models are shown with red dashed, green dotted, and magenta dot-dashed points. The 
horizontal dashed lines indicate the flux errorbars in 60-100 MeV band from the left panel. 


disk radii 0 < i? < Rmax with Rmax, which shrinks with 
the energy, reaching 0.15° in the 10-300 GeV energy bancQ 
see corresponding panel of Fig. The Rmax value is fixed 
by the requirement that the TS value decreases by 9. This 
source extension corresponds to the 3cr upper limit on the 
source size. 

The first six panels of Fig. also show the dependence 
of the estimate of the source flux in each energy bin on the 
assumptions about the spatial extent of the source. One 
could see that the flux estimates vary by no more than 
10% — 20% for the source extensions within Rmax- 

The bottom right panel of Fig |^ shows the 3a upper 
limits on the source extension as a function of energy. We 
find that the upper limits on the disk size become weaker 
with decreasing energy. This is obviously explained by the 
energy dependence of the PSF. The upper limit on the 
source size can be described with the single power law 
Rrnax — 0.44° (F/lGcV)with the precision ^ 15% at 
each point. 

In GeV range the fit prefers the diffuse source (of ^ 0.2° 
radius) t o the point source. This fi n ding is in the agreement 
with e . g. [Hooper fc Linden (|2011 |; Hooper & Goodenough 
(|2011 1; |Abazajian et al.| (|2014|), who hnd that the data is 
better described by a slightly extended source. However, 
evidence for the source extension is weak, below the 3a 
level. 


2.1.2. Fermi/LAT spectral analysis 

The spectral model of each source in the region of interest is 
chosen to be a power law with the slope —2 in each energy 
bin. The diffuse background model spectra are following the 
templates, the iso_source_v05_revl.txt for the isotropic 
background and the gll_iem_v05_revl. f its template for 
the Galactic diffuse emission background. The analysis 
was performed with the BinnedAnalysis and UpperLimits 
python modules provided with Fermi Science Tools soft¬ 


ware. The resulting spectrum of the source is shown in Fig. 

H 

In addition to statistical errors, we include a 10% sys¬ 
tematic err or in energy range 100 MeV - 1 TeV (Ackermann 
et al.||20T2|. This systematic error characterises uncertain¬ 


ties in the knowledge of the Fermi/LAT telescope. Taking 
the complexity of the GG region into account, we perform 
a separate study of additional systematic errors stemming 
from the uncertainties of the knowledge of properties of 
the interstellar medium (and as a consequence, of the dif¬ 
fuse Galactic emission) and the point source distribution in 
the region of interest. We do this by repeating the analysis 
using several templates for Galactic diffuse emission and 
several sets of point sources. For the Galactic diffuse back¬ 
ground, we consider gll_iem_v05_revl. f its, the most re¬ 
cent template, recommended Fermi/LAT collaboration for 
the analysis gll_iem_v05. fits , which is an earlier version 
of the template, and LRYusif ovXC04z6R30_Tsl50_mag^ 
a template based on the prediction from GALPROP 


code (Vladimirov et al. 2011). As a radical alternative to 
these models, we also consider a model in which we re¬ 
place the background in the 3° x 1° box around the GG in 
gll_iem_v05_revl. f its template with a rectangle of uni¬ 
form surface brightness. We consider the flux of this rect¬ 
angle to be an additional fitting parameter in each energy 
band. We verify that all the four models yield fluxes of the 
GG consistent within 10%. Ghanging the limiting flux level 
for the catalogue sources, or changing the catalogue alto¬ 
gether from 2FGL to 3FGL also produces changes in the 
flux that are within the 10% level; see right panel of Fig. 

The broadening of the Fermi/LAT PSF from ~ 0.1° at 
100 GeV to ^ 10° at 100 MeV together with the overall 
complexity of the region can lead to the source confusion 
problem. The spectral results in this case are given by the 
sum of the fluxes of several nearby point sources and/or the 
unaccounted diffuse emission. The variation of the number 
of model sources as well as the consideration of the different 


^ This size is similar to the distance (0.093°) between the GC 
and the closest to its source in the 3FGL catalogue. 


® See e.g. http://fermi.gsfc.nasa.gov/ssc/data/access/lat/ 
Model_details/FSSC_modeLdiffus _reprocessed_vl2.pdf 
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background models show that the flux attributed to the 
GC varies within 10% systematics in the 60 — 100 MeV 
band. This reveals that the possible presence of extended 
emission does not influence the measurement of the GC 
source spectrum. 

Energy dependence of the Fermi/IAT PSF could lead 
to yet another problem for the spectral extraction if the 
GC source is extended. The region from which the source 
flux is collected might have larger spatial extent at low 
energies because of wider PSF. To estimate the influence 
of this effect on the measurement of the spectrum, we have 
forced the spectral extraction for an extended source of the 
radius i? = 1°, which corresponds to Rmax at 100 MeV. 
This spectral extraction method ensures that the source 
flux is collected from the same region (rather than from 
within the PSF) at all energies. A comparison of the source 
spectrum extracted with this approach, with the spectrum 
obtained for the point source at the position of the GG, is 
shown in the bottom left panel of Fig. The two spectra 
are consistent with each other. This shows that the energy 
dependence of the PSF does not have a strong influence on 
the measurement of the spectrum. 


2.1.3. Fermi/LAT timing analysis 

The 1-10 GeV light curve, produced using the likelihood 
analysis, is shown in Fig The source does not exhibit 
variability higher than ~ 10% on yearlong timescales, the 
light curve is consistent with a constant flux model. 


2.2. INTEGRAL/ PICsIT data analysis 


We complement the GeV Fermi/LAT observations with 
data from MeV-range PIGsIT instrument. The Pixellated 
Imaging Gaesium Iodide Telescope (PIGsIT) is the high- 
energy detector of the IBIS telescope on board the 
INTEGRAL satellite, operating in energy range between 
175 keV and 10 MeV ( Di Cocco et al.|[2003 Labanti et al. 


20031. 


For the INTEGRAL/ PIGsIT analysis, we selected all 
data obtained from February 2003 to August 2012. The 
data consists of detector images obtained in eight energy 
bands for single and multiple events. The resulting effective 
exposure times on the GG are 8.6 and 12.4 Msec for sin¬ 
gle and multiple events, respectively. Detector uniformity 
and background maps have been built every ten spacecraft 
revolutions (three days each). 


The main sources of non-uniformities are two electronic 
effects relevant in the selection of single and multiple events 
on board. Multiple events are not detected at the bound¬ 
aries between semi-modules and for the pixels located out¬ 
side of the 4-by-4 pixel detector elements. The ratio be¬ 
tween the single and multiple events is therefore not uni¬ 
form across the detector. In addition, delays of time coin¬ 
cidence windows depend on individual front end electronic 
chips. Once corrected for these effects, the sky images are 
perfect but the instrument responses have to be corrected 
to take the average energy dependence of the multiple event 
selection into account. As this effect cannot be calibrated 
properly, the responses were adjusted using spectra of the 
Grab nebula. 


Mosaic images were finally built from the individual 
sky images. No source is detected at the Galactic centre 
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Fig. 4. The light curve of the Galactic centre in 1-10 GeV 
energy band. Each point represents the bin with duration of 
1 yr. The statistically allowed fits with constant and sloping 
lines are shown with dot-dashed curves. 
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and upper limits were derived from the variance images 
and converted in flux units assuming a power-law spectral 
model. Upper limits on the source flux derived from the 
INTEGRAL data are shown in Fig. 


3. Discussion 


The origin of the multiwavelength emission from the 
Galactic centre has been extensively discussed in the litera- 


& Neronov ( 

2005a); 

& I'akahara 

(Wi) 




see e.g 

. Aharonian 

(|2011 

I; 

Kusunose 


them can be divided into two broad classes, leptonic and 
hadronic. 

The spectrum of the GG shown in Fig. [3] is consis¬ 
tent w it h the previously reported spectrum of|Chernyakova| 


et al.| ( |2011| ), but extends to lower energies (down to 
60 Me V). 'I'his potentially opens a possibility of distinguish¬ 
ing between the leptonic and hadronic models because the 
hadronic spectral models are generically expected to pos¬ 
sess a spectral feature in the 100 MeV energy range (energy 
comparable to the mass of the pions). 

We obtained the spectrum of the GC under the assump¬ 
tion that the diffuse background in the region is reasonably 
well described in one of the templates above for the Galactic 
diffuse emission. These templates rely on the model of 
three-dimensional density distribution of the interstellar 
medium (except for the “constant rectangle” model tem¬ 
plate). The largest uncertainties of reconstruction of this 
three-dimensional distribution from different radial velocity 
data is in the direction of the Galactic centre. This uncer¬ 
tainty propagates to the uncertainty of the knowledge of the 
spatial morphology of the Galactic diffuse emission in the 
Galactic centre region and further to the uncertainty of the 
measurement of the flux of the Galactic centre source. In 
this respect, we are reassured by the fact that the measure¬ 
ment of the flux based on the “constant rectangle” model, 
which completely ignores the details of the spatial morphol¬ 
ogy of the diffuse emission around the source, agrees with 
the measurement based on the “standard” diffuse emission 
templates. 
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E,eV 

Fig. 5. Electron energy-loss cooling times for the medium 
parameters similar to those observed in the vicinity of the 
GC. The timescale of 300 yrs, corresponding to the flaring 
activity age, is shown with a thin solid line for convenience. 


The uncertainty of the three-dimensional matter distri¬ 
bution in the direction of the Galactic centre also leads to a 
potential problem of “source confusion”. The non-variable 
source in the direction of the Galactic centre might be con¬ 
fused with an unrelated molecular cloud, occasionally pro¬ 
jected to the GC position. In this case, a part of the ob¬ 
served source flux should be attributed to the “passive” 
molecular cloud source rather than to a source actively ac¬ 
celerating particles in the Galactic centre. The spectrum of 
this passive source should approximately repeat the spec¬ 
trum of diffuse emission from the interstellar medium (more 
precisely, of its pion decay and bremsstrahlung compo¬ 
nents). A test for this hypothesis could be a comparison 
of the spectral shapes of the Galactic centre with the inter¬ 
stellar medium around the Galactic centre. Fig. shows a 
comparison of the spectrum of the Galactic centre with the 
spectrum of the rectangular box of the size 3° x 1° (grey 
shaded region) from the “constant rectangle” diffuse back¬ 
ground model. The spectrum of the box is normalised to 
match the flux of the GG at ^ 300 GeV. One can see that 
the presented spectra are different. The spectrum of the 
GG in the energy band above 1 GeV is significantly softer 
than the spectrum of the diffuse emission right around the 
Galactic centre. This indicates that the flux of the source 
is perhaps not dominated by the flux of a passive molec¬ 
ular cloud positionally coincident with the Galactic centre 
on the sky. If a significant fraction of the source flux still 
comes from a molecular cloud situated in the GG region 
close or around Sgr A*, the mass of the cloud should be 
comparable to the entire mass of the nearby Sgr B2 com- 
plex( ^ IO^Mq ), which is also a source of high-energy 
and very high-energy, gamma-ray emission with the flux 
comp arable to that of the G alactic centre at 100 GeV en¬ 
ergy (Aharonian et al. 2006). Under this assumption, the 
upper limit on the size of the GeV source (< 0.1°) translates 
to a lower limit of > 10 ^ cm^ on the density of the medium 
within a 10 pc scale region. This is somewhat higher than 


A complex (see e.; 

y. Tsuboi et al. 1999 

Goldsmith et al. 

1990 

Genzel et al. 

2010 . 


Apart from the passive gamma-ray source produced by 
the interactions of interstellar cosmic rays with high-density 


medium near or around the Galactic centre, an active accel¬ 
eration process injects high-energy protons and/or electrons 
into the medium. Interactions of these particles lead to the 
gamma-ray emission of “leptonic” or “hadronic” origin. 


3.1. Leptonic model 


The MeV-TeV spectrum of the GG can be interpreted with 
a simple one-zone leptonic model if one takes the flaring na¬ 
ture of the source into account. In this case the emission is 
a combination of the emission from electrons injected dur- 
ing the strong flare, which occurred ~ 300 yrs ago (see e.g. 
Sunyaev et aL|[T9^ Koyama et al.||1996 2008 1 and during 


much weaker, recent activity. Electrons propagate through 
a medium with the density of about 10 ^ cm“^, through 

at 


,-3 


~ 0.3 — 3 eV (Davidson et al. 

1992 

Mezger et al. 

1996 

Hinton & Aharonian |2007 ) and magnetic held of ore 

er 10 - 


energy-loss channels are ionisation, bremsstrahlung, inverse 
Compton (IC), and synchrotron mechanisms. Fig. ^shows a 
comparison of the cooling times for these energy-loss chan¬ 
nels, as a function of electron energy. 

We find that the observed source spectrum could be 
well reproduced by models with rather different choice of 
parameters. Depending on the ambient medium density, the 
observed spectrum can be explained either with a pure IC 
model (low medium density, see left panel of Fig. m, or 
by a combination of a IC scattering with a bremsstrahlung 
emission. 


In the left panel of Fig. ^ we show an inverse Compton 
dominated model fit in which electrons are scattered on 
photons with energy 0.5 eV and energy density 5 x 
10'^ eV/cm^. The magnetic held is assumed to be 2.5 x 
10“'* G. Electrons are injected with a spectral index -1.5 for 
10 yrs during the hare that occurred 300 yrs ago. After the 
hare, the luminosity decreases by factor of 1.7 x 10^. During 
the hare the GC has reached a luminosity of ~ 2 -10^® erg/s 
in the 7 -rays (if the electron spectrum during the hare ex¬ 
tended up to the 10 TeV energy band). 

The observed GeV bump in the GC spectrum is ex¬ 
plained by IC scattering of the soft photons on the electrons 
injected during the hare. The cooling time of electrons with 
energies about 10 GeV is comparable to the time since the 
hare. Suppression of the hux above the GeV energy is ex¬ 
plained by the fact that electrons emitting at higher ener¬ 
gies have already cooled down. The emission above 10 GeV 
is mostly due to the electrons injected during the ongoing 
low activity period of the source. Above 10 TeV the emis¬ 
sion is suppressed by synchrotron losses, which exceed the 
1C losses in the Klein-Nishina regime of 1C scattering. To 
explain the TeV spectral shape in y-rays, the magnetic held 
should be - 10“'* G. 


The inverse Compton model of the GeV bump of the 
spectrum can be verihed through the variability properties 
of the signal. If the GeV bump is really due to the electrons 
that were injected 300 years ago during a major hare of 
the source, gradual cooling of these electrons should lead 
to the displacement of the bump towards lower energies in 
the future. A noticeable displacement should already occur 
on the 10-yr timescale. The 1 — 10 GeV hux is expected to 
decrease by ~ 5 — 10% on this timescale. From the source 
light curve shown in Fig.|^one could see that the sensitivity 
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Fig. 6 . Left panel: X-ray to TeV energies sp ectrum of the Galact ic centre. The data are adopted from Baganoff et al. 
(|2003|) (green solid cross, X-ray quiescent), |Porquet et al.| (|2008|) (green dashed cross. X-ray flaring), this work (red 
points) and Aharonian et al. ( 2009|) (blue points, HESS). The iJN'flEGRAL/ PICsIT upper limits (this work) are shown 
in black. The leptonic model proposed in this work for the low-density case is shown with the solid black line. Right 
panel: the black line shows the hadronic model of 7 -ray radiation from the broken power-law distribution of protons. 
The model significantly under predicts the flux at < 150 MeV energies; see text for details. The statistical error is shown 
with thick errorbars, while the thin ones indicate systematic levels, which are shown in the left panel. 


of Fermi/LAT is marginally sufficient for the detection of 
this decrease. 


3.2. Hadronic model 


Aharonian & Neronov (2005b|) and jChernyakova et al. 
( |201ip developed a model in which the GC Fermi/HESS 


spectrum is explained by hadronic radiation from relativis¬ 
tic protons, which diffuse by ^ 5 pc away from the cen¬ 
tral source. Data above 100 MeV are well described by 
the model for different sets of parameters (density profile 
of surrounding medium, characteristics of the diffusion co¬ 
efficient, injection rate history). A generic feature of the 
hadronic models is a low-energy cut-off in the spectrum at 
100 MeV at pion production threshold; see right panel of 
FigUl 

In this figure, we show the results of the fitting of 
Fermi/LAT data with the model incorporating broken 
power-law proton spectrum at Fermi energies , used as an 
approximation of Chernyakova et al. (20111 model. We 
find that this model provides a reasonably good fit of the 
data {x^/NDF = 11.0 /22, with included systematics, and 
/NDF = 32.2/22, without). The two lowest-energy data 
points (60—100 MeV and 100 — 150 MeV), where the cut-off 
in hadronic models is expected, deviate from the model by 
2 . 4(7 and 1.4 ct (3.5(t and 2.2 ct, without systematic errors). 
In absolute values, the model flux in these energy bins is 
^ 5 and ~ 1.7 times lower than the observed flux. 

The observed discrepancy with the data can be ex¬ 
plained by the emission of several leptonic mechanisms, 
possibly operating in the region. The bremsstrahlung or 
IC emission from the electrons in the region can perfectly 
fill the gap in this energy interval. 

4. Conclusions 

We have presented the results of ~ 6 years of Fermi/LAT 
observations of the Galactic centre. We have shown that 
the spatial morphology of the GC source is consistent with 


the point-like source. The 3 (t upper limit on the radius of 
the source was found to be ~ 0.13° at energies > 10 GeV, 
and ~ 0.22° at energies 3—10 GeV. For lower energies, the 
radius of the source does not exceed ~ 0.7°; see Fig. 

The spectrum of the source shown in the left panel of 
Fig. [^covers a broader energy range co mpared to the previ¬ 
ous work by Chernyakova et al. (2011). In the overlapping 
range (0.1-3(10 GeV), these two spectra are consistent with 
each other. This indicates that the spectrum is indepen¬ 
dent of the version of the Fermi/LAT data analysis software 
as well as of the choice of the diffuse/point source back¬ 
grounds, confirming the correctness and robustness of pre¬ 
viously obtained results. The spectrum in the 60-300 MeV 
energy band does not show any evidence for low-energy 
cut-off, which is expected in pure hadronic models of the 
Galactic centre emission. 

The MeV-TeV band spectrum of the GC is well mod¬ 
elled with a leptonic model in which the GeV bump in the 
spectrum is produced by the IC emission from electrons 
injected during the strong flare from the GC, which hap¬ 
pened ~ 300 yrs ago. The model has a testable prediction of 
a decrease of the 1-10 GeV flux by 5 — 10% on the ~ 10 yrs 
timescale. The quality of the Fermi/LAT data is insufficient 
for testing the model prediction with a 6 -yr data set (see 
Fig[^. The prediction, however, can be tested within the 
next decade before the end of Fermi mission, or w ith one 
of its succe ssors, e.g ASTROGAAj^ (d AMMA-40 0 (iGalper 


et al. |2014D. HERD.(|Zhang et al.pOldD, PANGU (|Wu et al. 
2014P , or DAMPEfj 


Also, our analysis does not rule out the models in which 
the leptonic GeV component of t he GC spectrum has a 
differ ent origin fro m the leptonic (Kusunose & Takahara 
2012p or hadronic (|Guo et ^|2013P TeV component. 


The hadronic models (Aharonian & Neronov 2005b 
Chernyakova et al. 2011) conflict with the data. The ac¬ 


count of the electron bremsstrahlung or IC emission can 


® http://astrogam.iaps.inaf.it/Doc/extract.pdf 
^ http://dpnc.unige.ch/dampe/index.html 
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increase the flux below ~ 100 MeV and removes this con¬ 
flict. We conclude that the presence of a leptonic component 
in the cosmic rays near the GC is unavoidably required in 
order to match the predictions of the models with the data. 
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